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Starting from a microscopic description of a system of strongly interacting electrons in a strong 
magnetic field in a finite geometry, we construct the boundary low energy effective theory for a 
fractional quantum Hall droplet taking into account the effects of a smooth edge. The effective 
theory obtained is the standard chiral boson theory (chiral Luttinger theory) with an additional 
self-interacting term which is induced by the boundary. As an example of the consequences of this 
\ model, we show that such modification leads to a non-universal reduction in the tunnelling exponent 

which is independent of the filling fraction. This is in qualitative agreement with experiments, that 
systematically found exponents smaller than those predicted by the ordinary chiral Luttinger liquid 
theory. 

as : 

■ PACS numbers: 



INTRODUCTION 



Quantum Hall effect (QHE) has been the subject of intense research since 1980 when the first plateaux in the Hall 
conductivity were observed for integer filling fractions The interest on these systems was then reinforced after 
the fractional QHE (FQHE) was identified as a novel state of matter 0. Since then, different approaches have been 
developed and applied to the study of these systems, such as the successful variational ground-state wave function 
proposed by Laughlin 3] . This description led to the concept of fractionally charged quasiparticles as a consequence 
of the strong correlations in the system, which were in turn experimentally observed [J]. The composite fermion 
approach [j| , which provides a description of the FQHE as a QHE where a quasiparticle is built up as a bound state 
of an electron and a given number of magnetic flux quanta, also provides a consistent description. An alternative 
description of these systems, based on the Chern-Simons topological action also reproduces many of the experimentally 
observed features M. 



in 

oo 

In spite of the enormous success of all these descriptions, a theory for the QHE derived systematically from a 
microscopic description is still lacking. It is the purpose of the present paper to give one step in this direction and 
provide an effective description of the FQHE which is derived from a first principles microscopic action describing 
interacting electrons in a strong external magnetic field in a controlled way. In doing this, we are able to complement 
previous descriptions and provide a way to take into account the effects of a smooth boundary. These effects have 
been already considered in 0, EH EL EH EH EH EH , the main difference with our approach is that our starting point 
. £h ! i s the microscopic model and hence it does not rely on a phenomenological description. 

The study of tunnelling of electrons into incompressible and compressible quantum Hall states has been the subject 
of intense research from both the theoretical and experimental sides (see Q and references therein) . It is found that 
the tunnelling conductivity is non-Ohmic, / oc V a , with a being a function of the filling fraction v. This behavior 
can be understood within the chiral Luttinger liquid description of the edge physics advanced by Wen 9] . Although 
tunnelling experiments in fractional quantum Hall effect (FQHE) systems have shown certain degree of agreement 
with the theoretical predictions obtained from the chiral Luttinger liquid picture @ , there remain certain discrepancies 
which have been addressed by different authors @, EH EH EH EH EH EH- In particular, a reduction of the order 
of 10% from the theoretical prediction for the tunnelling exponent has been observed experimentally, as well as the 
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absence of the theoretically predicted plateaux structure of this exponent (Tol , [T3L [l5j . These issues have been the 
subject of intense debate Q- 

Here we present an alternative derivation of the boundary effective action, by including higher order terms in the 
low energy expansion which are induced by the smoothness of the edge. This paper can be considered as a formal 
derivation of the model first presented in [8[. In our calculations, the presence of a smooth boundary shows to be 
crucial e. g. in the computation of tunnelling exponents. The effects of the edge shape and the confining potential on 
these observables have been recently addressed using numerical exact diagonalization fl2l |. where it was found that 
the result could deviate from the chiral Luttinger liquid universal predictions. Our results are in agreement with these 
findings. 

In the simplest case of an almost sharp edge, the effective action we obtain corresponds to a self-interacting chiral 
boson. Using this improved effective action we compute the tunnelling exponent which, to one loop order, receives a 
small non-universal negative correction depending on the electron density and the short distance cutoff (the "effective 
size" of the particles). 

The paper is organized as follows: In Section 2, we present the model which corresponds to a microscopic description 
of a fully polarized self-interacting electron gas in a transverse magnetic field. The action is invariant under relabeling 
of particles, a symmetry that plays an important role in the continuum limit (corresponding to the low energy 
description) as we show in Section 4. In Section 3 we rewrite the action by making a boost transformation to 
eliminate the external electric field and we discard terms quadratic in the particle velocities. This last approximation 
is justified due to the presence of a strong external magnetic field. We show that as a result of the interactions, the 
system is projected into the fundamental state and that this projection corresponds to the imposition of a constraint 
fSection l3.2p which defines the effective degrees of freedom as the fluctuations along the flat directions of the potential. 
In Section 4 we construct the low energy effective action by going to the continuum limit and show that the resulting 
action is topological. We show that the original symmetry under relabeling of particles corresponds to the invariance 
under area preserving diffeomorphisms of the effective action. A condition which arises from the solution of the 
constraint is that our analysis is applicable to incompressible states (Section |4~4"1) . Although the constituent particles 
are fermions, the collective degrees of freedom, which correspond to fluctuations along the flat directions of the 
interaction potential, are described by bosonic fields. 

In Section 5 we first solve the constraint by performing a gradient expansion on the coordinate fields and then, 
by boosting back to the original (laboratory) coordinate system, we arrive at one of our main results, which is the 
effective action (|5.5I) . In Section 6 we look in more detail into the effective action close to a smooth boundary, which 
requires to recover the discrete structure along the direction perpendicular to the boundary, as expected. This leads 
generically to a boundary effective theory of AT coupled chiral bosons, with M = W/a, a being the mean inter- 
particle distance and W the effective width of the edge (|6.1[) . By integrating out the internal degrees of freedom, 
we finally obtain the effective action for the boundary modes, which corresponds to that of a chiral self-interacting 
bosonic field. From this action, we compute the correction to the tunnelling exponent, the result being in qualitative 
agreement with experiments. The result obtained for the exponent is non-universal, which is a direct consequence 
of the inclusion of a smooth edge [l6l |. A higher loop computation together with a renormalization group analysis 
should be done, but this is out of the scope of the present paper. The main aim here is to obtain for the first time 
an effective field theory, starting from a microscopic description of the real problem and which includes the effects of 
smooth boundary, a feature present in almost all experiments. The observed edge reconstruction effect [lij could be 
in principle recovered within our approach by including the effects of the fluctuations fields perpendicular to the flat 
directions of the interaction potential (field f2 in (|4.12[) h This is the subject of future investigations. We conclude in 
Section 7 with a discussion of the results and open perspectives. 

2. POLARIZED ELECTRON GAS IN A TRANSVERSE MAGNETIC FIELD 

2.1. The action 

The most general action for a system of N interacting electrons is given by 

S = J dt {^j4 2 ~ V KK]j * = 1 ' 2 P,9=1,-",N , (2.1) 

where the degrees of freedom of the system are represented by the 2N functions x£(f), i.e. the positions of the electrons 
in a cartesian coordinate system. These functions will be our dynamical variables. For our purposes it is enough to 
consider spinless fermions in two space dimensions since, as explained later, typical experiments are performed in the 
presence of strong magnetic fields. 
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Note that this action is completely generic and makes no simplification. The electrons interact with the background 
and with each other by means of a general potential U[x p ,x p ]. This potential may depend simultaneously on the 
positions and velocities of all particles (i.e. all p = 1, • • • ,N), and to keep this in mind we write its argument in 
square brackets. 

We split the interaction into an external part and an electron-electron interaction part U[x p ,x p ] = V ext [x p ,x p ] + 
V mt [x p ,Xp]. Since electrons interact independently with the external field, the external part can be written as a sum 
of a single term per electron V ext [x p ,x p ] — J2 P V ext (x p ,x p ) and we then have 

s = J dt {j2 (y4 2 - ^(4>4)) - v mt KK^j > ( 2 - 2 ) 

where the first term corresponds to the kinetic contribution. 

Since the system is non-relativistic, the external potential V ext (x p ,x p ) can be expanded in powers of velocities as 

^(4,4) = u(4) + eA,(4)ij + \5m{x> p )x^ + 0{x^) (2.3) 

the first two terms representing the electric and magnetic coupling to the external field. The quadratic order in 
velocities represents a position dependent correction to the mass that can be reabsorbed in the kinetic term, while 
higher orders have been discarded. 

On the other hand, the expansion of the inter-particle interaction U mt [x p ,Xp] reads 

N r)V int N N Fp\fint 

v mt 4,4] = v mt 4, o] + J2 -^j- 4- °H + E E 4= °i H x * + o& 3 ) (2.4) 

The resulting terms can be grouped according to the power of velocities involved, resulting in a mass like (particle 
mixing) term, a magnetic term linear in velocities, and an electric potential term 



f)2trint \ N ( f)V int \ N ^ 



(2.5) 



This is the action of our system. Up to this point no approximation has been made (other than the non-relativistic 
limit) and the action has been kept completely general. 



2.2. An important symmetry 

Since the particles are identical, the action is naturally invariant under arbitrary permutations of the particle 
indices. We will call this property invariance under relabeling of the particles. This can be formally formulated as 
follows: the redefinition 

p - p' = P '(p) (2.6) 

acts in the dynamical variables as 

4 = x p(p') = X V (2-7) 

leaving the action unchanged. 



3. THE INTERACTIONS AND THE RESULTING PROJECTION INTO THE FUNDAMENTAL STATE 

The general expansion for the external potential 12.3|) resulted in an electromagnetic interaction in three space-time 
dimensions. We define the corresponding external magnetic field as B = Bz, with B = diAje^ and we assume it is 
constant and homogeneous in all space to conform with the typical experimental setup. Moreover, we take it to be 
very strong so as to completely polarize the particle spin in the z direction, justifying in this way our choice of spinlcss 
fermions to represent the electrons. The electric field on the other hand is defined as Ei = —diU. It is assumed to be 
zero in the interior of the region of space occupied by the sample, and smoothly growing at the edge of that region. 
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3.1. A convenient coordinate system 

In this setup, external electric and magnetic fields are perpendicular and, in the absence of additional interactions, 
we can discard the acceleration term in the equations of motion = Ei + e^v^B ~ 0. Solving for we see that 
the resulting classical particle motion takes place in the direction normal E^. With an appropriate boost, we could 
jump to a coordinate system in which the particle is at rest and there is no electric field. Even in the presence of an 
inter-particle interaction as in our case, we are interested in eliminating the external electric field by such a boost. 
The corresponding velocity will be 

Vi = ~€ijdjU (3.1) 

Using the invariant expression B 2 (l — (diu) 2 /B 2 ), and expanding in powers of E/B, we can solve for the magnetic 
field in the boosted coordinate system B as 
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The action in the new coordinate system reads 



p I 1 N I ppxrint \ u N ( r)V int \ \ 

S = dt \ ?£ \(m + 6m(x p ))S pq — — S^"E ^(5$) + — — x 3 p -V mt (3.3) 

J \ M V dx p dx qJ P V dx p I J 

where Aj(St l p ) is the external vector potential in the new coordinates. Note that in this frame u(x p ) = and the 
external magnetic field B(x p ) is non-homogeneous. In other words, the effects of the external electric field in the old 
coordinate system x, is now encoded in B(x p ) in the new system x p . This is very important since, as we will see, it 
turns out to lead to localization of the excitations. 

On the other hand, it is evident that the relabeling symmetry persists in the new variables, taking the same form 
as before, i.e. 

K = K' ( 3 - 4 ) 



3.2. The projection constraint and a new form of the action 

Due to the presence of a strong external magnetic field, we can discard the kinetic terms quadratic in velocities, as 
well as the magnetic inter- particle interaction dV mt /dx p . In a perturbative quantization approach, this approximation 
corresponds to the projection into the lowest Landau level (LLL). The scrupulous reader can keep all the terms, then 
go through a quantization procedure and project the result into the LLL. 

The resulting effective Lagrangian in then linear in velocities and the resulting action for the degrees of freedom in 
the ground state is 



Hi 



N 

£ e i,(x;)x J p -u mt [ x ;]] , (3.5) 



f [ X p] = V* nt [x£,0]. To be consistent, we should supplement this action with the condition 
enforcing the low energy projection, namely that the system remains in its ground state 7i = E m i n . Within our 
approximation, this can be written as 

V lnt [x p ]=Q (3.6) 

where we have eliminated any non-vanishing contribution to the ground state energy by a shift in the potential. This 
constraint can be enforced at the action level by introducing a Lagrange multiplier X(t) to obtain 

s = idt\y ]eA 3 (x;)x j p xv mt [x;^j , (3.7) 

the equations of motion being now obtained by varying the action with respect to x p and the Lagrange multiplier 
X(t). This new form of the action is still invariant under the relabeling symmetry. 
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4. THE CONTINUUM LIMIT AND THE RESULTING TOPOLOGICAL FIELD THEORY 

We assume that all the important physical scales involved in our problem are large compared to the inter-particle 
distance. These scales can be characterized in terms of the relevant length of the external field B/djB. Since the 
derivatives involved in this definition can differ significatively when taken in different directions, so does the accuracy 
of this assumption. In the present section we assume that it is valid in any direction, an assumption which is well 
justified in the bulk. We will come back to this point in the case in which this assumption ceases to be valid, i.e. 
when analyzing boundary effects in Section [6] 

4.1. The continuum space 

We take the continuum limit in the above system as follows: first we define our continuum space and functions 
living on it by the following steps: 

1. We introduce an auxiliary two-dimensional space parameterized by variables y\ 

2. We set a (not necessarily regular) lattice y p in that space by the condition y p = x^(0), that will be useful later. 

3. We use the discrete quantities x p (t) of our problem to define lattice functions x l (y p , t) on the auxiliary y p space 
by means of 

x*(g,f) = (4.1) 
In other words, we label each particle p with their initial position y p . 

4. We re-interpret these lattice quantities as the values at the lattice points of some continuum fields x 8 (jp',t) 

x i (g,t)=x i (yi,t)| y<=4 (4.2) 

These continuum fields can be understood as an interpolation of the discrete quantities representing the particle 
positions at a given time in the intersticial region of y-space. There are many interpolation methods that could 
be used for that [13], however, we do not need an explicit prescription for it. Note that this re- interpretation is 
consistent with the initial condition on the continuum variables 

x*(j?,f = 0)|- J= -, = x»(y£,0) = x*(0) = f p (4.3) 

4.2. The resulting field theory 

We reformulate now the dynamics in terms of the continuum variables introduced in (|4.1[) by starting from the 
discrete action (|3.7p . The procedure goes as follows: 

1. We make a simplicial decomposition of the y l plane consistent with the lattice, i.e. having a single point y p for 
each simplex. It is easy to convince oneself that this is always possible and that in two dimensions the simplices 
are triangles that can be defined as having a lattice point at each vertex. The easiest choice is a simplicial 
decomposition topologically equivalent to the triangular lattice. 

2. We call AA p the (small) area of the p-th simplex and multiply and divide by it the p-th. term in the action (|3.7|) . 
This results in a Riemann approximation for the integral 

N N 1 

J2eMm,t))dtm,t) = ^Al p ^ e I i (x(y;,f))a t x(g,f) = 

p p P 

~ J d 2 yp(y)eA 3 (x(y\t))d t x(y>,t) + 

where the number density function p(y) stands for the interpolation of p(y p ) = 1/ AA p . This approximation of 
the sum remains accurate as long as the jump on the argument when going from one lattice point to another 
stays small. We have kept terms up to linear order in these jumps and provided an estimation of the error made 
in that process that will be useful later. 



di lepA^dtx) Ay p _ AA P (4.4) 
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Applying the above steps explicitly we obtain the continuum low energy effective theory 

S = Jdt(J d 2 ypeA J {x l )d t x- XV mt [x}^j , (4.5) 

in terms of the continuum fields x l (j/,i) and the Lagrange multiplier X(t). 

We are interested in solutions x l (y,t) of the field theory (|4.5[) that are good approximations of the true solutions 
x p (t) of the discrete microscopic system. In particular this means that the error estimated in (|4.4p has to be small 
when evaluated on them. In other words, we should discard those solutions whose derivatives are large compared to 
the inter-particle distance. 

Note that the action contains a local velocity-dependent term, coming from the external magnetic field, plus a 
non-local potential constraint, coming from the inter-particle interaction. The topological nature of the resulting field 
theory is evident in view of the absence of any spatial derivatives. 

In this part of the calculation a new function p(y) has appeared. We need some suitable ansatz for the form of this 
function. Since is represents the particle number density in the ground state, a judicious ansatz would be to choose 
it consistently with the forces on these particles in that state. In the y plane, particles at rest experience no forces 
other than their repulsive interaction, so we can safely choose p(y) as a constant p Q in this plane. 



4.3. Symmetry under Area Preserving Diffeomorphisms and gauge invariance 

The relabeling symmetry (|2.6p - (|2.7|) . (13. 4[) is inherited by the continuum limit theory, where it manifests as a 
symmetry under area preserving diffeomorphisms. To see this, we first note that the area AA p of the simplex 
associated to the particle p becomes, after a relabeling, 

Ai p = Ai p(p0 = Ai;, (4.6) 

The second equality is evident from the fact that AA P becomes, after relabeling, the area of the simplex associated 
to the particle p' . 

In the continuum limit, the symmetry can be written as 

y ^ y' = y(y) (4.7) 

while the fields transform according to 

x i '(y',t)=x i (y,t) (4.8) 

The equality (|4.6p then implies that 

d 2 y = d 2 y' (4.9) 

Then the transformation corresponding to the relabeling of particles becomes an area preserving diffeomorphism in 
the continuum theory 

Due to the definition of p, the combination pd 2 y is invariant under any transformation. Then the kinetic term in 
the action is obviously invariant under area preserving diffeomorphisms. The potential, on the other hand, inherits 
its invariance from the discrete case. 

We conclude that our system has a continuos invariance. In the case in which B is constant, this corresponds 
to the gauge invariance of the usual Chern-Simons description [9J. On the other hand, in the case in which B is 
space dependent, the only transformations leaving the action invariant are those which are time independent. This 
implies that gauge invariance is broken in our theory, in a similar way as the imposition of a boundary breaks gauge 
invariance in [9j, enforcing the introduction of the boundary chiral boson. In our formulation this chiral boson theory 
arises naturally as a consequence of the low energy constraint. This is the first insight of one of the main claims of 
this paper, i.e. that a smooth boundary, represented by an external potential varying in space, can be described by 
a theory whose qualitative behavior is similar of that of the chiral boson. 



4.4. Flat directions of the ground state and a simpler form of the constraint 



The Lagrange multiplier in the action imposes the constraint that forces the system to be in the ground state. To 
obtain the remaining degrees of freedom we would then need to solve this constraint. This seems impossible without 
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knowing explicitly the form of the potential, but we show in this section that this is indeed simplified by the symmetry 
under area preserving diffcomorphisms of our system. In consequence we can replace the constraint for an explicit 
one, without giving any details on the potential. 

The condition (|3 .€>[) fixes all directions except those that leave the potential invariant. In other words, particles 
can move along the flat directions of the potential, while staying in the ground state. These directions correspond 
then to the remaining degrees of freedom. To determine them, let us note that any small deformation x l (y,t) — 
x* (jj, t) + Sx l (y, t) around a given minimum of the potential x l (y, t), can be written in the form 

Sx l (y,t) = M tj (y,t) (fyfifo t) + e jk d k A(y, t)) (4.10) 

for suitable (fi, A), when M is a given invertible matrix field that we choose for convenience as My = diX 3 (y, t). With 
this choice we can write 

5x\ Vl t) = d l xi(y,t)(d j n(y,t) + e jk d k A(y,t)) (4.11) 

or in other words 

?(y,t) = x'fe- + fyfifot) + e jk d k A,t) (4.12) 

implying that any small deformation of the configuration x 3 (y,t) can be obtained by a small diffeomorphism on y. 
The f2 part represents dilatations of the y plane, while the A part corresponds to area preserving deformations. Going 
now to the potential, this implies 

V[x {y, t) + 5x?{y, t)} = V [x (y, + d,Q + e jk d k A, t)} 

where in the last equality we use the invariance of the potential under the area preserving diffeomorphism generated 
by A, 8%ji = eijdjA. We conclude that the condition (|3.6| fixes and leaves A completely free. The flat directions of 
the potential are then the variations of x l (y,t) generated by area preserving diffeomorphisms. 

In particular, since x l (y, 0) = y 1 , we conclude that the time evolution map, going from y % to x 2 (y,t), is area 
preserving. In other words, the Jacobian of the time evolution is unity, implying that a physical configuration satisfies 

|^'| = 1 (4.14) 

This constraint has a simple and intuitive interpretation: note that the particle density (in x l space) at time t is given 
in terms of the initial density p Q by the formula 

**.*)= |^| (4-15) 

i.e. our evolution is that of an incompressible fluid. For this reason, we will concentrate in what follows in in- 
compressible states. To describe compressible states on the other hand, we would have to consider the f2 degree of 
freedom. 

It should be clear now that we can replace the constraint (|3.6p by the more explicit one (|4.14[) . In other words, we 
can replace (I4.5[) by the completely local action 

S = Jdt J d 2 y~ Po (e A^dtx-X (|c^|-l)) (4.16) 

where now X(y,t) is a local Lagrange multiplier imposing (|4.14|) . 

As a conclusion, in this section we have been able to use the symmetries to describe the low energy dynamics 
independently of the explicit form of the inter-particle potential. 

4.5. The conserved charge and an immediate consequence 

Since the action is invariant under the group of area preserving diffeomorphisms in the y l plane, we can calculate 
the corresponding Noether charge. To do that, we write the infinitesimal transformation as 



y'i=Vi + ^jdjA 



(4.17) 
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and the charge reads 

Q=\ J #V dj&d&K (4.18) 

Moreover, since this charge is conserved for any function A, its integrand is conserved, i.e. the magnitude 

B(x)e ab e lj d j x b d l x a = 2B{x)\d i x i \ = 2£(x) (4.19) 

is conserved. Here in the last equality we used the constraint (|4. 14j) . We can fix its value according to the initial 
condition as 

B(x) = B(x) = B{y) (4.20) 
t=o 

Note that this relation implies that the value of the magnetic field at the particle position is the same along all its 
motion. In other words, the particle moves along the level lines of the function B. Again, this may be interpreted as 
a behavior analogous to that of the Chern-Simons theory, in which excitations move along the boundary, and can be 
described by a chiral boson. Since in our case the boundary has been replaced by a smooth space-dependent external 
field, this translates into a motion along the level lines of this external field, and the corresponding chiral boson theory 
will necessarily have to take this into account. 

5. SOLUTION OF THE CONSTRAINT AND THE RESULTING BOUNDARY THEORY 
5.1. Solution of the constraint as a gradient expansion 

The next step is to solve the constraint (|4. 14|) and rewrite the action in terms of the dynamical degrees of freedom. 
Here and in what follows, we use as parameters the filling fraction in the x plane, defined as v — p/eB, and the 
effective (dimensionfull) parameter 9 = l/(2irp ) defined in ,18| as the square of the non-commutativity length. 

If we parameterize the x 1 field with the help of a slowly varying bosonic field ip according to the following gradient 
expansion 

x k = f + 8e kl d l cb+ ( ^ e kl e ij d^d^ + O(0 3 ) (5.1) 

then (|4. 14(1 is solved up to order 8 3 . This field represents the effective degrees of freedom after projection into the 
ground state. It can be interpreted to first order in 9 as representing the coordinate fluctuations. Incidentally we see 
in (|5.ip that <f> is defined up to the addition of an arbitrary function of time cj> —></> + g(t), which is the usual gauge 
invariance of the chiral boson theory. 

Replacing this solution in the action (|4.16[) we get 

S = ^2 Jdtd 2 yi>-h l ^d ] (0d l d t (P) + p ] (e kl d k dt(f ) d l <j ) dic^ + 0{6 2 ) (5.2) 

5.2. Boost back 

Note that the variables x l (y,t) refer to the position of the particle initially located at y. This is the Lagrangian 
description of a fluid, where each fluid clement is labeled with its initial position and then followed along its motion 
throughout the plane. But the initial position y l is given in an awkward coordinate system that is related to that 
of the laboratory by the transformation discussed in Section 13.11 We would like to refer things to the laboratory 
coordinate system by going back on the original frame as y l = y r (y). Note that for convenience we do not apply 
this operation to the x l fields. They still correspond to the position at time t in the tilted system, but now the fluid 
element is identified by its initial position y l on the laboratory frame. 

Under this transformation, the action becomes 

J dtdPyv- 1 c« (d^diidt - v a d a )cj>) + ^d j (e kl d k (d t - v a d a )<t> dicj>)\ (5.3) 
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where we defined the position dependent filling fraction in the (physical) y plane as v — p/eB, and it satisfies 

Po( 1+ 4) 



eB V 2B 2 

On the other hand, the velocity of the excitations is defined as in eq. (|3.1[) 

In the special case in which B is constant the action is a total derivative. The usual procedure is to assume this 
and then impose a boundary to the region in which the fluid moves. We choose a different way here, allowing the 
parameter v to change in space and introducing through it the information about the geometry, without adding any 
boundary. A further integration by parts gives 



S = J dtd 2 yf(^cl > d l (d t -v a d a )cj ) +^e bc d b (d t -v a d a )cbd t cl > d c 



(5.5) 



where we defined the vector field U(y), tangent to the level line of v at y, as 

f = jtdsiu- 1 ) (5.6) 

We see that the degrees of freedom are localized at the positions at which the change of variables discussed in Section 
13. H is non-trivial, i.e. where the external electric field diu is non-constant. Moreover they propagate in the direction 
of ti along which E 2 is constant. 

Here we see that, wherever the sample is homogeneous, the derivatives in front of each term vanish, leaving us 
without any dynamics. On the other hand, the dynamical degrees of freedom localize at places where there is a 
change in the properties of the material (i.e. the filling fraction). This is precisely what happens at the edge of the 
sample, and it is at the core of our argument. 

To stress this point, let us suppose that the filling fraction changes as a step function, then their derivatives will 
provide delta functions leading us to a one dimensional boundary theory, which corresponds to a Chiral Luttinger 
Liquid description [^] . Our treatment however is more general in the sense that it includes the possibility of dealing 
with smooth edges. 

The action (|5.5p is written completely in terms of the dynamical degrees of freedom <fr. Note that this is not a 
boundary theory because it is defined in the full two dimensional space. Nevertheless, the degrees of freedom are 
bounded to the region in which there is a change in the parameter v and, as predicted, propagate chirally along their 
level surfaces. 

6. RECOVERING THE DISCRETE STRUCTURE NEAR THE EDGE 

When we approach the edge, the derivatives of the functions entering into the action are no longer negligible. Then 
the assumption we made when we took the continuum limit (|4.4p is not satisfied, and the solutions of the continuum 
theory do not represent faithfully the underlying discrete system. To go around this problem, we need to recover the 
discrete structure near the edge. Calling yu the direction that runs along the edge and y± the perpendicular direction, 
we see that the y± derivatives become big in the edge region. Since y\\ derivatives are still small, we can keep the yu 
integral while in the direction y± transverse to the edge, we proceed as follows: 

1. Cut the space in the y± direction in N slices at y\ of width Aj/j_ ~ a. 

2. Replace the integral by a sum over the values of the integrand evaluated at each slice, and define 4> n (y\\) = 

<t>(y\\,vl)- 

3. Replace all the y± derivatives by its finite difference approximation 

dxttvw , y±) = (<t>(y\\ , vl) - 4>{y\\ , vT l ))l*v±. 

With this method we obtain an effective theory for the fields <f> n (yu) living on each slice n, with action 

S = -g-2 j dt d V\\ ^ ( K n^r((/>n) + «n (AsQftn) + £r(<ftr» <ftn-l)) ) (6-1) 

71 J 71 = ^ ' 
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where we have defined the free, self interaction and interaction Lagrangians as 

£f{4>ti) = <t>nd\\{d t -V n d\\)(j) n 



Cs{(t>n) = 4>n{d t - V n d\\){d\\(/> n ) 2 
£x (</>„, 0„-i) = -<t> n -x(d t - v n d\\)(d\\4>„) 2 (6.2) 



and the constants K n and k' u are given by 

11 ,6>/l 



(6.3) 



It is important to note that the field <j) n will enter into the action only when the constants n ni n' n are non- vanishing, 
i.e. if there is a change in the properties of the sample between the slices n and n + 1. 
The gauge invariance is now <j) n — ^ cf> n + g(t) adding the same git) to all the 0's. 

Note that when the sample has a sharp edge, i.e. if the density changes within a region whose width W is smaller 
that the slicing length Ay±, i.e. W < Ay± ~ a, then the whole procedure is not applicable and the boundary theory 
corresponds to the usual chiral boson theory Q. 

Let us suppose that the edge is wider than the slicing length, so that we have to keep a finite number of terms of 
the sum in (|6.ip . n = 1, • • • ,N. 

We have then a single non-dynamical field ip = <P-i i n the interaction term £j(0 o , </>-i), all other fields appearing in 
the Cx{4>m 0n-i) terms being dynamical. The integration of this multiplier field will enforce a constraint on the field 
<Pq to which it is coupled. When solved, 4>o = f{y\\ + v$t), and replaced in the action, it implies that the interaction 
term with </>i takes the form 

- ^<fdtd y] \f(y\\ +v t)(d t -«id||)(c>||<M 2 (6-4) 

which is different from zero provided that the two velocities are not equal. Since we have assumed that the electric 
field in the edge zone is approximately constant in regions of size a ~ Ayx, hence vq w v\ and we can discard this 
term within the present approximation. 

We have then obtained an action very similar to the original one but without any Lagrange multiplier 



1 / . . . v-^ 



8tt 2 



n=l 



—r / dtdyu y K„£^(^„) + K n (C S ((j>n) + (^nl - l)^x(0n>^n-l)) (6-5) 



where the lower bound of the sum has changed and the coefficient (S n i — 1) ensures that there is no interaction term 
for the field <\>\ . 

7. AN APPLICATION OF THE MODEL: TUNNELING EXPONENTS 

Let us consider the case of an almost sharp boundary, in which the density can be well approximated by choosing 
J\f = 1. In that case the action is 

5 =-^2 Jdtdy^K^d^dt-vd^ + n'^idt-vd^d^f) (7.1) 

where we have dropped the subindexes since we have a single field. This is a chiral boson theory with a cubic higher 
derivative self interaction, the latter arising essentially from the smoothness of the edge. To see the effect of this last 
term in the physical properties, we compute the propagator since it is directly related to the tunnelling exponent 0]. 
In Fourier space we write for the Feynman propagator 

G(p, Up ) = — -— (7.2) 

Kp (Ul p — Vp) + It 

while the vertex is proportional to 

- k'5(p + q + r)(uj p - vp) qr (7.3) 



11 



Then, the one loop corrected propagator reads 



2ir 

G'(P,U P ) = —, r—^f r (7.4) 

up [Up - vp) - dG(p,u) p ) 



We consider those Feynmann diagrams whose contribution to the corrected propagator is 



k' 2 

SG(p, lj p ) = 7rA 2 — (u) p - vp)p (7.5) 



where A is a momentum cutoff. Replaced in (|7.2p this gives 



G(p,u>p) = - — 1 (7.6) 

[k + dn)p [Up — vp) 



where the non-universal correction to the level is given by 



k' 2 

5k = -7T— A 2 (7.7) 



Note that the effective coupling is k' ' / ' K 2 n = 9 2 /4Ax\, being independent of the details of the boundary. A non 
vanishing correction to the velocity arises from the remaining diagrams but it produces no physical consequences. 

Since the tunnelling exponent is directly related to the level by a « n, the appearance of a correction for the latter 
implies, to first order, a correction to the tunneling exponent. It should be stressed at this point that the correction 
is negative, in qualitative accordance with the experimental results. This is the main result of this section. 

To estimate the magnitude of this correction, we need to relate the momentum cutoff A to the minimal space 
distance measurable a, which naturally leads to the identification A = it /a. Using (|6.3[) we then get for this choice of 
the cutoff 



(4a z p y 16 

where we have further identified p a = 1/a 2 . The predicted dependence on the density of this non- universal correction 
could in principle be tested experimentally. 

Putting it all together this leads to a linear dependence of the tunneling exponent a as a function of \jv which 
is very close to the experimental fit presented in [10( (see for example Fig. 3 of this reference). However it should 
be kept in mind that our present result is valid for incompressible situations. It can be observed that the departure 
of our result from the experimental data becomes greater for increasing l/i/, which could be attributed to the need 
to work with TV > 1 in (|6.5p . This could in turn be related to the need to consider a wider boundary region. The 
analysis of the higher M case will be presented separately . 



8. CONCLUSIONS 



To summarize, we have constructed an effective field theory that describes the boundary region of a fractional 
quantum Hall effect droplet, starting from the microscopic description. We included the interplay between the 
granularity of the fluid (the inter-particle distance a) and the effects of a smooth boundary. This extended boundary 
has been defined by means of a non- vanishing electric field, that leads to a space-dependent filling fraction within this 
region. The physical excitations of the system are localized at the boundary, and are described by a system of self 
interacting chiral bosons. 

All the approximations we made along our calculations are under control: by taking the non-relativistic limit in 
the strong magnetic field regime, we enforced the projection into the LLL. Keeping the interaction between particles 
completely generic, this low energy projection results in a constraint on the dynamical degrees of freedom of the 
model. With the help of the symmetries of the system, we were able to replace this constraint by an equivalent one 
that implies incompressibility and is completely independent of the details of the interaction potential. 

As an example of the possible applications of our approach, we have computed the correction to one loop order of 
the tunnelling exponent to lowest order in the boundary width, which is in qualitative agreement with experimental 
results. 
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The results presented here can in principle be improved in a systematic way. Wider boundary regions can be easily 
treated within the present approach, by simply taking higher values of J\f in (|6.5[) . Higher loop computations could 
also be envisaged, however they would involve the solution of the constraint to higher orders in 0. This is out of 
the scope of the present paper. As a further improvement, one can describe compressible situations by including the 
effects of the degrees of freedom transverse to the flat directions of the potential (described by f2 in (|4.13p ). 

Other approaches to include a boundary in a Chern-Simons description of a quantum Hall droplet have been 
proposed [20(. In this context, it would be interesting to study their connection to our model. 

It should be mentioned that the experimental data in [Io| was reanalyzed in fl3j in view of the results of fl5j , and 
certain degree of agreement with the plateau structure was obtained. Here we propose an alternative description 
which provides a reasonable agreement with the raw experimental data as presented in [Toj . It would be interesting to 
study the higher order theory that we have constructed in the present paper along the lines of the approach presented 
in ref. [l5l ]. This could provide a closer agreement between theoretical and experimental results. 

It would be also interesting to analyze the consequences of the higher order corrections induced by the smoothness 
of the boundary that we have obtained here in the transition between plateau states. 
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